Phase separation and crystallisation of polydisperse hard spheres 
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Hard spheres with a polydispersity above approximately 8% are shown to crystallise into two 
phase-separated solid phases. A polydispersity above 8% is too large to be tolerated by a single solid 
phase but phase separation produces two fractions with polydispersities sufficiently narrow to allow 
them to crystallise. It may not be possible to observe this in experiment due to the intervention of 
a glass transition. 



PACS: 82.70.Dd, 81.30.-t, 64.70.Dv 

Experiments on synthetic colloids are inevitably on 
polydisperse colloids. The colloidal particles are not all 
the same, some are larger than average, some are smaller 
101 . This has prompted considerable theoretical study of 
the phase behaviour of polydisperse systems, in particu- 
lar the fluid-solid transition of polydisperse hard spheres 
|^|-^. The conclusion of these studies is that a solid of 
hard spheres can tolerate a polydispersity of a little less 
than 0.1, beyond this it is never stable with respect to the 
fluid phase. The 0.1 refers to the ratio of the standard 
deviation to the mean, denoted below by s. As Pusey 
has noted, this makes sense if we note that the lat- 
tice constant of the solid is at most 10% larger than the 
diameter of the hard spheres. Hard spheres melt at a 
fraction ~ 0.74 of close packing ||^ which corresponds to 
a lattice constant ~ 1.11 times the hard sphere diame- 
ter. Only if the polydispersity is sufficiently narrow that 
few or no spheres are larger than this lattice spacing can 
the mixture crystallise 0. We have derived an analyt- 
ical free energy for the solid phase of polydisperse hard 
spheres and our results agree with those of previous au- 
thors (2|H] . We find that the Helmholtz free energy of the 
solid phase of polydisperse hard spheres with s > 0.082 
is never lower than that of the fluid phase at the same 
density. 

However, all this takes no account of possible demix- 
ing. If the spheres phase separate into fractions then 
these fractions will have narrower polydispersities than 
the original mixture and so can crystallise. Of course, 
this phase separation costs some ideal mixing entropy but 
at high density spheres pack more efficiently in the solid 
phase than in the fluid resulting in a larger entropy in the 
solid. The possibility of demixing has been mentioned be- 
fore but no calculations have been performed. Here, 



calculations show polydisperse hard spheres crystallising 
into two solid phases, one containing the larger spheres 
and one containing the smaller spheres. There is no fluid- 
fluid phase separation at the moderate polydispersities 
we consider, although there may well be at much larger 
polydispersities ||^,|lo|]. By moderate we mean polydis- 
persities up to about 20%. 

The crystallisation of polydisperse hard-sphere-like 
colloids into more than one solid phase has not been ob- 
served experimentally. This is despite the fact that it is 
almost certainly the equilibrium behaviour. The discrep- 
ancy between theory and experiment may be due to the 
intervention of the glass transition [ pTp^ ]. The cost in 
demixing entropy and the less efficient packing of polydis- 
perse spheres into a solid vis-a-vis monodisperse spheres 
both act to raise the coexisting densities of the fluid-solid 
transition. It is quite possible that when the polydisper- 
sity is sufficiently large for the equilibrium behaviour to 
be crystallisation into two solid phases, the transition 
occurs at densities above the glass transition and is thus 
impossible to observe in experiment. Already for near- 
hard-sphere colloids with a polydispersity s = 0.075 the 
fluid-solid and glass transitions appear to be very close 

0. 

In a polydisperse mixture of spheres, spheres with a 
range of diameters are present ||l^,|l^; there is a con- 
tinuous distribution of sizes of spheres with a number 
density px{a)d(T of spheres of size a. p is the total num- 
ber density of spheres. The width of the polydispersity 
is characterised by a width parameter w. The larger w is 
the broader the distribution of sizes present in the mix- 
ture. In the limit w ^ we recover a monodisperse 
system. The functional form of the polydispersity should 
not matter too much; we select a simple form, the hat 
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function: 








(7 < — w/2 

a - w/2 < <T <a + w/2 
a > 7j + w/2 



(1) 



where a is the mean diameter. The standard devia- 
tion divided by its mean, of the distribution x{a) is 
s = w/{VT2a). 

An approximation to the Helmholtz free energy of the 
sohd phase may be obtained using a ceU theory |15[| . Our 
system is athermal so the free energy is equal to minus the 
entropy (ignoring the contribution from the momenta). 
Within a cell theory the entropy per particle is just the 
logarithm of the configuration space available to one par- 
ticle trapped in a cell formed by its neighbours. These 
neighbouring particles are taken to be fixed at their lat- 
tice positions, which are the lattice constant a from the 
centre of the cell. For a solid of monodisperse spheres of 
diameter a the hard-core repulsion prevents the particle 
from moving more than (a — a) from its lattice position 
in the direction of one of its neighbours. Then we can 
approximate the cell by a sphere of radius a — a [l5| ] and 
so obtain a free energy per particle fs in the solid phase 



fs 1 



47^/ ^3' 
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where for a close-packed lattice, a — [rjcp/ilY^^] k is 
Boltzmann's constant and T is the temperature. The 
volume fraction rj of monodisperse spheres is related to 
the number density by 77 = {Tr/6)pa^, and its maximum 
value is rjcp = 7r/(3\/2), which is achieved at close pack- 
ing. 

In order to generalise Eq. (||) to a solid of polydisperse 
spheres we have to make a further approximation. We 
neglect any correlations between the size of spheres at 
adjacent lattice points. In reality we would expect that 
on average smaller spheres will tend to cluster around 
large spheres to reduce the strain on the lattice. Every 
sphere is assumed to be in a cell formed from spheres 
of the average diameter a. Then the cell of a sphere of 
size a is of radius (a — (cr + ct)/2). This gives us the free 
energy of a sphere of size a; we simply integrate over all 
sizes of spheres, weighting with x{<t) to obtain the free 
energy per particle 



kT 



= -ln(47r/3)-3 / x{(j)\n[a- (a + ct)/2] da - Inw 
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where the last term is the ideal mixing entropy p3[ for 

The lattice constant 



our polydisperse mixture |13 



a is still derived from a — {rjcp/rjY^^ with rjcp for a 



monodisperse solid. For polydisperse spheres with a dis- 
tribution Eq. (]l|) the volume fraction is related to the 
total number density by 77 = (7r/6)p(l 4- w^/4). For the 
free energy of the fluid phase we use the polydisperse gen- 
eralisation of Salacuse and Stell of the free energy of 
Boublik, and Mansoori, Carnahan, Starling and Leland 
(BMCSL) 

Finding phase coexistence for polydisperse mixtures is 
very involved but see Refs. [[8|jl|]. Determin- 

ing whether the fluid or the solid phase is more stable, 
has the lower free energy, at a given density is however 
straightforward, so we will content ourselves with this. In 
Fig. 1 we have plotted the free energies of the fluid and 
sohd phases with a polydispersity w/a = 0.3 (s = 0.087). 
The free energy of the solid phase is always higher than 
that of the fluid phase and so is never the equilibrium 
phase. This holds if and only if the polydispersity is 
greater than w = 0.28 (s > 0.082). However, the free en- 
ergy of a solid phase with half the width, here w/a = 0.15 
does drop below that of the fluid phase. As the fluid 
has a higher free energy than solid phases with half the 
polydispersity it is unstable with respect to them. The 
phase separation will not of course be that complete or 
simple. Describing the two phases as one with all the 
spheres smaller than a and the other with the spheres 
larger than ct is a gross oversimplification. However, we 
have proven that the fluid becomes unstable with respect 
to two solid phases, whilst being stable with respect to 
one. We have also plotted the free energy of the solid 
phase with w/a = O.l. In the range of densities plotted 
in Fig. 1 it is always higher than for w/a = 0.15 so we 
expect that on crystallisation phase separation into two 
not three solid phases will occur. Also, the free energy of 
the fluid phase with w/a — 0.3 is lower than that of the 
fluid phase for either w/a — 0.15 or 0.1. As an aside, we 
note that polydisperse sticky spheres also phase separate 
into more than one solid phase 20 1. 

The free energy of phase separated solid phases only 
drops below that of the fluid phase above a volume frac- 
tion rj = 0.6. Experiments on near-hard-sphere colloids 
with a polydispersity s w 0.04 [|l2| show a glass tran- 
sition at 7/ — 0.58. Although of course we have not 
been able to find the coexisting fluid and solid phases 
this at least suggests that, in experiment, the combined 
phase-separation-crystallisation phase transition may be 
prevented by a glass transition. Even if it is not, the 
dynamics of a combined phase separation and crystalli- 
sation will be slow 0j2^ . 
^^•j Hard spheres with a polydispersity width w > 0.28 
(s > 0.082) cannot crystallise into a single solid phase, 
its free energy is always higher than that of the fluid, 
(^feowever, they can crystallise into two solid phases, we 
have verified this for polydispersities up to w = 0.4 
(s = 0.12). We have not performed calculations beyond 
this because our solid phase free energy may not be re- 
liable at these polydispersities, and because the density 
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at which the fluid and solid free energies are equal are so 
high (77 > 0.63) at these polydispersities that the BMCSL 
free energy of the fluid phase may no longer be reliable. 

In conclusion, experiments on near-hard-sphere col- 
loids only find crystallisation when the polydispersity is 
no more than approximately s — 0.075 This is not 

the equilibrium behaviour of polydisperse hard spheres. 
At equilibrium, as the polydispersity is increased a sec- 
ond solid phase appears; we predict that this occurs at 
a polydispersity s < 0.082. Each of the two solid phases 
is about half as polydisperse as the original fluid. For 
polydispersities s of tens of per cent the picture is much 
less clear: there will now be spheres with very different 
diameters present which may favour the formation of lat- 
tice types other than a substitutionally-disordered close- 
packed lattice. The formation of very different lattices 
has been observed for binary mixtures of hard spheres 
|^,|3|. Another possibility is that highly polydisperse 
hard spheres do not have an equilibrium solid phase. 
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Fig. 1. The free energy of polydisperse hard spheres 
in the fluid and solid phase for several polydispersities. 
The solid curve is for the fluid phase with a polydisper- 
sity width w — 0.3. The dotted, dashed and dot-dashed 
curves are for the solid phase with polydispersity widths 
w ~ 0.3, 0.15 and 0.1, respectively. The dotted curve has 
been cut off just before a rapid divergence. 
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